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Abstract 

Solar proton events (SPEs) are known to have caused changes in constituents in the Earth’s polar neutral middle 
atmosphere. The past four years, 2000-2003, have been replete with SPEs and huge fluxes of high energy protons 
occurred in July and November 2000, September and November 2001, and October 2003. The highly energetic 
protons produce ionizations, excitations, dissociations, and dissociative ionizations of the background constituents, 
which lead to the production of HO x (H, OH, HQz) and NO y (N, NO, NO 2 , NO 3 , N 2 0 5 , HNO 3 , HO 2 NO 2 , CIONO 2 , 
B 1 -ONO 2 )- The HO x increases lead to short-lived ozone decreases in the polar mesosphere and upper stratosphere 
due to the short lifetimes of the HO x constituents. Large mesospheric ozone depletions (>70%) due to the HO x 
enhancements ware observed and modeled as a result of the very large July 2000 SPE. The NO, increases lead to 
long-lived stratospheric ozone changes because of fee long lifetime of the NO, family in this region. Polar total 
ozone depletions >1% were simulated in both hemispheres for extended periods of time (several months) as a result 
of fee NO, enhancements due to the very large SPEs. 


1. Introduction 


Explosions on the Sun sometimes result in large fluxes of high-energy solar protons at the Earth, especially 
near Solar Maximum. This period of time, wherein the solar proton flux is generally elevated for a few days, is 
known as a solar proton event (SPE). Solar cycle 23 experienced a large number of extremely energetic SPEs in 
years 2000-2003. Huge fluxes of high-energy protons occurred in July and November 2000, September and 
November 200 1 , and October 2003 . 

Solar protons are guided by the Earth’s magnetic field and impact both the northern and southern polar cap 
regions (>60° geomagnetic latitude), e.g., see Jackman and McPeters (2004). These protons can impact the neutral 
middle atmosphere (stratosphere and mesosphere) and produce ionizations, dissociative ionizations, and excitations. 
Both HO x (H, OH, H0 2 ) and NO, (N, NO, N0 2 , N0 3 , N 2 0 5 , HN0 3 , H0 2 N0 2 , C10N0 2 , BrON0 2 ) constituents are 
produced either directly or through a photochemical sequence (e.g., S wider and Keneshea, 1973; Crutzen et al. 1975; 
Jackman et al. 1980; Solomon et al. 1981; McPeters 1 986; Zadorozhny' et al. 1992). Ozone is also impacted by the 
solar protons through direct photochemical destruction forced by the HO x andNO y enhancements (e.g.. Weeks et al. 
1972; Heath etaL 1977; Solomon etal. 1983; Jackman etal. 1990). 

The SPEs that occurred in years 2000-2003 were noteworthy and details of their atmospheric response, 
including both satellite measurements and model predictions are included in this discussion. The paper is divided 
into six primary sections, including the introduction. We discuss the very important solar proton measurements and 
their production cf odd hydrogen (HO x ) and odd nitrogen (NO,-) in section 2. A comparison of the SPEs in solar 
cycle 23 with some of the largest in past solar cycles is also undertaken in section 2. The GSFC two-dimensional 
model used to simulate the impact of the SPEs on the atmosphere is discussed in section 3. The short-term impact of 



these SPEs on ozone during and for several days after particular events is given in section 4. Longer tom influences 
of the SPEs on the middle atmosphere are discussed in section 5. Finally, the conclusions are given in section 6. 


2. Proton Flaxes; Odd Hydrogen (HO,) and Odd Nitrogen (NO y ) Production 

Solar proton fluxes are measured by a few satellites in interplanetary space or in orbit around the Earth. The 
most accessible and useful proton flux data are available from the National Oceanic and Atmospheric Administration 
(NOAA) Space Environment Center (SEC) for the NOAA Geostationary Operational Environmental Satellites 
(GOES) [see http:// sec. imaa.gov/Data/goes.html] . GOES proton fluxes are provided in several energy intervals (>1 
MeV, >5 MeV, > 10 MeV, >30 MeV, >50 MeV, and >100 MeV) at this site, updated every five minutes. GOES-8 
data are considered most reliable for proton fluxes depositing energy into polar latitudes (private communication, 
Terry Onsager, NOAA SEC). GOES-8 data are, therefore, used for the periods January 1, 2000 to April 8, 2003; 
and May 10, 2003 to June 18, 2003. GOES-1 1 became the primary satellite for protons on June 19, 2003 and was 
used as the proton flux source through December 31, 2003. GOES-10 data was used to fill in the gap of missing 
proton flux data from April 9 to May 9, 2003. 

The solar proton fluxes were used to compute daily average ion pair production profiles using the energy 
deposition methodology discussed in Vitt and Jackman (1996). Odd hydrogen (HO,) is formed through complicated 
ion chemistry (Solomon et al. 1981). Each ion pair is assumed to produce two HO, constituents up to an altitude of 
approximately 70 km. Above 70 km, the HO, production is assumed to be that provided by Solomon et al. ( 1 981 , 
Figure 2). The HO, constituents have lifetimes of only hours in the middle atmosphere, therefore, any further effects 
on other constituents from the HO, group are apparent only during and shortly after an SPE. 

Atomic nitrogen is produced by the primary protons and associated secondary electrons causing 
dissociations, predissociations, or dissociative ionizations in collisions with N 2 . Following Porta - et aL (1976) and 
Jackman et aL (1980), we assume that 1.25 N atoms are produced per ion pair. The N atoms rapidly produce NO 
and other odd nitrogen (NO y ) constituents. Odd nitrogen has a relatively short lifetime (-days) in the sunlit middle 
and upper mesosphere, however, lower mesospheric and stratospheric NO y can last for weeks past an SPE. A mostly 
dark middle atmosphere in the late fall and winter conserves a large portion of the SPE-produced NO y , which can 
then be transported to lower altitudes via the general downward flowing winds during this time of year. The lifetime 
of this enhanced NO y can range from months to years if transported to the middle and lower stratosphere. 

We have quantified middle atmospheric NO y production before (Jackman et al. 1980, 1990; Vitt and 
Jackman, 1996) for years 1955 through 1993. We add NO y computations in this study to these earlier calculations 
for years 1994 through 2003 and present the annual production from SPEs for the 49-year period 1 955 through 2003 
in Figure 1 . The source of proton flux data for years 2000 through 2003 was explained above. For the years 1 994 
through 1999 we use two satellites: 1) GOES-7 for the period January 1, 1994 through February 28, 1995; and 2) 
GOES-8 for the period March 1, 1995 through the end of 1999. The annual-averaged sunspot number is also shrwn 
in Figure 1 to illustrate the rough correlation between solar maximum periods and frequency of SPEs. 

Solar cycle 23 was quite active with several very large SPEs (also, see Krivolutsky et al. 2003), especially in 
years 2000 (July and November), 2001 (September and two in November), and 2003 (October). Substantial 
amounts of NO y were produced in 2000, 2001 , and 2003 with 1989 being the only year showing a larger production. 
The ten largest SPEs in the past 40 years are given in Table 1 with six of them occurring in the past solar maximum 
period. 


3. GSFC Two-Dimensional Model 

The latest version of the Goddard Space Flight Center (GSFC) two-dimensional (2D) atmospheric model 
was used to predict atmospheric changes caused by the solar protons. The model has been in use since the late 
1980’s and has undergone extensive improvements over the years (Douglass et al. 1989; Jackman et al. 1990). 
Fleming et al. (2002) describes the methodology to compute the transport for the GSFC 2D model. This techniqic 
uses the global winds and temperatures from the United Kingdom Meteorological Office (UKMO) data assimilatbn 
system for the years 1992-2000. The photochemical reaction rates have been recently updated to Sander et al 
(2003). The GSFC 2D chemistry solver has been improved and now uses the Atmospheric Environmental Reseach 
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(AER) 2D model scheme (Ko et aL 1984). The new chemistry solver allows a diurnal cycle simulation every day 
and provides for a more accurate simulation of atmospheric constituents. 


4. Short-term Atmospheric Influences from Solar Proton Events 

We used the GSFC 2D model to compute two primary simulations, the “base” and the “perturbed” 
simulations, for the years 1998 through 2005. The transport for years 1998-2000 is driven by the UKMO products 
for those particular years, whereas the transport for the individual years 2001-2005 is a repeat of that derived for 
2000. The “base” simulation includes no SPEs, whereas the “perturbed” simulation includes all SPEs from January 
1, 2000 through December 3 1, 2003. The perturbation to the atmosphere was caused by the SPE-produced HO* and 
NOy enhancement 

We show the HO* production (in # cm'V 1 ), ozone depletion (in percent), and NO* (N0+N0 2 ) enhancement 
(in ppbv) in Figure 2 as a result of the extremely large solar proton event that started on Bastille Day (July 1 4) in 
2000. The model computations are compared with Upper Atmosphere Research Satellite (UARS) HALogen 
Occulation Experiment (HALOE) ozone and NO* measurements. The HALOE observed constituent changes during 
July 14-18 were calculated by comparing to the background atmospheric amounts before the SPE, defined as the 
average of the July 12-13 measurements. The produced HO* constituents drive practically all the ozone depletion in 
the mesosphere and the upper stratosphere during the event There are several HO* catalytic destruction cycles for 
ozone. An example of one that is important in the middle and upper mesosphere is: 

H + 0 3 ^OH + Oz 
OH + O H +0 2 


Net: 0 3 + O Oz + Oz 

The NO* enhancement can cause some of the depletion in the upper stratosphere and lower mesosphere 
during the event, but forces all of the ozone depletion after July 16 (Figure 2c, 2d, and 2e). Ozone is reduced by the 
NO* constituents through the following primary catalytic destruction cycle: 

no + o 3 -»no 2 + o 2 
no 2 + o -» no + o 2 


Net: 0 3 + O 0 2 + 0 2 . 

The ozone changes are s imil ar between HALOE measurements and model computations during the 
maximum intensity of the July 14-16, 2000 SPE. Both measurements and model calculations show a maximum of 
about 70% depletion in the middle and upper mesosphere on July 15. The model computed ozone depletions in the 
lower mesosphere and upper stratosphere are larger than that measured by HALOE on July 16 (late in the day) and 
on July 17-18 (see Figure 2b and 2c). The ozone reduction by SPE-enhanced HO* is thus fairly reasonably 
simulated, whereas the ozone reduction by SPE-enhanced NO* appears to have difficulties. It is possible that the 
HALOE measured ozone change is underestimated due to some long-term secular change that may occur ova- a few 
days. 

Clearly, the HALOE and modeled NO* are enhanced in the mesosphere and upper stratosphere on July 17-18 
(Figure 2d and 2e) by file extremely large SPE. These NO* enhancements should lead to a reduced ozone, which is 
simulated by the model in Figure 2c. The computed enhanced NO* is somewhat larger than measured in the middle 
and upper mesosphere on July 17-18, but is similar to the observed NO* in the lower mesosphere and upper 
stratosphere during these two days. 


3 


5. Long-term Atmospheric Influences from Solar Proton Events 

Over the course of a few days, the enhancements of NO and N0 2 from SPEs elevate the amount of other 
NO y constituents. As noted above, the NO y family has a lifetime of months to years if transported to the middle and 
lower stratosphere. Transport to lower altitudes is especially effective in the middle to late Fall and Winter, when the 
middle atmospheric winds are directed polewards and downwards. 

Four very large SPEs occurred in the northern middle to late Fall time period (see November 2000, 
November 2001, and October 2003 in Table 1). The SPE-related NO y enhancements were thus produced at an 
opportune time to be conserved for a very long period of time. We computed the percentage change of NO y and 
ozone in the northern polar latitudes (50-90°N) for years 2000-2005 and present the results in Figure 3. NO y 
enhancements of greater than 100% are noted on several occasions in the upper stratosphere and lower mesosphere, 
however, the SPE-caused NO y increases in the middle to late Fall periods lead to much larger and longer-lasting 
middle and lower stratospheric enhancements. Much of the northern polar middle stratosphere has computed NO y 
enhancements of >10% in 2001, 2002, and 2004 as a result of the middle to late Fall very large SPEs. These periods 
are highlighted in gray in Figure 3 (upper plot). The increased NOy led to a northern polar stratospheric ozone 
depletion for extended periods. SPE-caused depletions in ozone greater than 2% are highlighted in gray in Figure 3 
(lower plot) to show the scope of this effect 

The impact on total ozone is shown in Figure 4. Total ozone is reduced by a maximum of about 3% in late 
2001 and 2002. Much of this depletion was caused by the huge amount of NO y produced by the July 2000 SPE in 
the Winter, which was conserved and transported downwards (Randall et al. 2001), where it steadily reduced ozone 
over a period of a few months. This July 2000 SPE-enhanced stratospheric NO y was slowly reduced over time, 
however, other large SPEs also added to the NO y reservoir in both the southern and northern hemispheres. Both 
polar hemispheres had extended periods of depleted ozone greater than 1% (hig h l igh ted in gray) as a result of SPEs. 


6. Conclusions 

Six very large SPEs occurred in solar cycle 23 and caused very si gnifi cant perturbations in the polar middle 
atmospheric regions. These SPE-caused influences were quite substantial (e.g., >70% ozone destruction), but also 
very short-lived (-days) in the middle and upper mesosphere. The SPE-produced HO* likely caused the ozone 
depletions in this atmospheric region. The NO y enhancements from the SPEs caused some ozone depletion on a 
short-time scale (days) during the events in the lower mesosphere and upper stratosphere. However, much of the 
NO y lasted for a longer time period of months past the very large SPEs. A longer-lived ozone destruction was 
connected to this SPE-enhanced NO y . Depletions greater than 1 % were computed for both northern and southern 
polar hemispheres. 
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Table 1 . Largest ten solar proton events in past forty years 


Date of SPEs 

Rankin 

size 

NO y production 
in the middle atmosphere 
(# of molecules) 

October 19-27, 1989 

1 

6.7 xlO 33 

August 2-10, 1972 

2 

3.6 x 10 33 

July 14-10, 2000 

3 

3.5 x 10 33 ’ 

October 28-31, 2003 

4 

3.4 x 10 33 

November 5-7, 2001 

5 

3.2 x 10 33 

November 9-1 1 , 2000 

6 

2.3 x 10 33 

September 24-30, 2001 

7 

2.0 x 10 33 

August 13-26, 1989 

8 

1.8 x 10 33 

November 23-25, 2001 

9 

1.7 x 10 33 

September 2-7, 1966 

10 

1.2 x 10 33 

*The number reported here for the July 2000 SPE is slightly 
larger than that reported in Jackman et al. (2001). GOES-8 proton 
fluxes were used in this work, whereas GOES-10 proton fluxes 
were used in Jackman et al. (2001). 
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Figure 1. Total number of NO y molecules produced per year in the 
polar stratosphere and mesosphere by SPEs (solid histogram - left 
ordinate) and annu a ll y-averaged sunspot number (dashed line - right 
ordinate) for years 1955 through 2003. 
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Figure 2. Polar Northern Hemisphere (near 65°N) 
pressure versus time cross sections during the perturbed 
period (July 14-18, 2000) for a) HO* production for 
contour levels 200, 500, 1000, 2000, 5000, 10000, and 
20000 cmrY 1 ; b) HALOE O 3 and c) model O 3 
decreases, both for contour levels -5, —10, -30, -50, 
and -70%; d) HALOE NO* and e) model NO* 
increases, both for contour levels 2, 5, 10, 20, 50, 100, 
and 200 ppbv. The HALOE ozone and NO* changes 
were computed by comparing to the background average 
of the July 12-13 observations. The model changes 
were computed by comparing the “perturbed” to the 
“base” simulation. 
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Figures. Model computed 
percentage changes in NO y and 0 3 
for the polar Northern Hemisphere 
area (50-90°N) for 2000-2005 
resulting from SPEs in 2000-2003. 
Contour levels for NO y (top plot) 
are 0, +1, +2, +10, +20, and 
+100%. The gray highlighted areas 
for NO y indicate increases greater 
than 10%. Contour levels for 0 3 
(bottom plot) are -20, -10, —2, —1, 
-0.2, 0, +0.2%. The gray 
highlight ed areas for 0 3 indicate 
decreases greater than 2%. These 
changes were computed by 
comparing the “perturbed” to the 
“base” simulation. 
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Figure 4. Model computed percentage 
total ozone changes from 2000-2005 
resulting from SPEs in 2000-2003. 
Contour intervals are -3, -2, -1 , —0.5, - 
0.2, and -0.1%. The gray highl ighted 
areas indicate total ozone decreases 
greater than 1%. These changes were 
computed by comparing the “perturbed” 
to the “base” simulation. 
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Brief, Popular Summary of the Paper: 

Several very large solar proton events occurred in the past solar cycle during years 2000- 
2003. These events caused changes in polar mesospheric and stratospheric ozone and other 
constituents. The influence of these natural processes on ozone need to be more folly understood 
so that the human-made impact on ozone can be more reliably predicted. 

Large-scale explosions on the Sun near solar maximum lead to emissions of charged 
particles (mainly protons and electrons) from the Sun. Occasionally the Earth is positioned in its 
orbit such that these solar particles interact with the Earth’s magnetosphere and rain down on the 
polar regions. “Solar proton events” have been used to describe these phenomena since the 
protons associated with these solar events cause the most atmospheric disturbance. 

Solar proton events create hydrogen- ami nitrogen-containing compounds, which can lead 
to ozone destruction in the mesosphere and upper stratosphere. Six of the largest ten events in 
the past 40 years occurred in the years 2000-2003. The solar protons connected with the July 
2000 event caused short-lived (-1-2 days) mesospheric ozone decreases of over 70% because of 
the increases in the hydrogen-containing constituents. Increases in the longer-lived nitrogen 
oxides caused several extended periods (-months) of polar total ozone depletions of greater than 
1 %. 



